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Abstract. The kinetics of the chiral phase transition is studied within a linear quark-meson- 
<r model, using a Monte-Carlo approach to semiclassical particle-field dynamics. The meson 
fields are described on the mean-field level and quarks and antiquarks as ensembles of test 
particles. Collisions between quarks and antiquarks as well as the qq annihilation to a mesons 
and the decay of a mesons is treated, using the corresponding transition-matrix elements from 
the underlying quantum field theory, obeying strictly the rule of detailed balance and energy- 
momentum conservation. The approach allows to study fluctuations without making ad hoc 
assumptions concerning the statistical nature of the random process as necessary in Langevin- 
Fokker-Planck frameworks. 


1. Introduction 

One of the motivations for the study of ultrarelativistic heavy-ion collisions is to gain a detailed 
understanding of the phase diagram of strongly interacting matter 1 . At the largest energies as 
achieved at the Large Hadron Collider (LHC) and the Relativistic Heavy Ion Collider (RHIC) 
a hot and dense fireball is formed which can be described to a surprising accuracy as a nearly 
perfect fluid of strongly coupled quarks and gluons (QGP) undergoing a transition to a hot 
hadron-resonance gas. In these situations, where the net-baryon density or the baryon-chemical 
potential are small, lattice-QCD (1QCD) calculations indicate a crossover transition from 
confined to deconfined matter as well as from a phase where chiral symmetry is spontaneously 
broken to one where it is restored at a (common) transition temperature T c ~ 160 MeV |2|. 

At lower collision energies, as studied in the RHIC beam-energy scan (BES) program, at the 
CERN SPS, and the future FAIR and NIC A experiments, the produced medium starts at lower 
temperatures and larger net-baryon densities. Since in this situation the application of 1QCD is 
challenging due to the “sign problem” at finite /xb, one relies on effective chiral models, which 
predict the existence of a first-order chiral-phase-transition line ending in a critical point of a 
second-order phase transition. 

For theory the challenge is to provide possible observables for this phase structure in heavy-ion 
collisions, like the (“grand-canonical”) fluctuations of conserved charges like net-baryon number 
or electric charge. Not only the question, how to effectively model the phase transition (e.g., 
with the Nambu-Jona-Lasinio (NJL) or the (linear) <r model with extensions taking into account 


gluonic degrees of freedom implementing Polyakov loops) arises but also, which of the features of 
the phase structure predicted for such models applying thermal quantum field theory (describing 
a medium in thermal and chemical equilibrium) like (critical) fluctuations of conserved charges 
survive for a rapidly expanding and cooling fireball as created in heavy-ion collisions. 

To address the latter question, one relies on transport simulations to describe the off- 
equilibrium dynamics of the fireball. One approach is the use of ideal or viscous hydrodynamics 
to describe the bulk evolution of the fireball (assuming a state close to local thermal equilibrium), 
which successfully describes key phenomena of heavy-ion collisions, and adding the fluctuations 
by hand in a Langevin approach [3j-|5|. On the other hand this implies that the statistics of 
the random process has to be put in as an ad hoc assumption. Usually a Gaussian Markovian 
(“white noise”) is assumed, but the simulation of non-Markovian (“colored noise”) processes is 
feasible in principle |6|. 

On the other hand, one would like to study the implication of different transition scenarios 
(cross-over, 1 st order, 2 nd order) as realized in the various quantum-field theoretical models 


on the nature of the hopefully observable fluctuations [7-13 on the statistics of the probably 
observable fluctuations. 

In this work we present a novel Monte-Carlo approach to address this challenging problem 
using the most simple quark-meson linear a model 114]. The meson fields are treated on the 
mean-field level and the quarks and antiquarks are realized in terms of a test-particle ensemble. 
Here the challenge is to implement “discrete” local interaction processes like elastic collisions and 
reactions like the qq annihilation to a a meson and the decay of a mesons to a qq pair admitting 
not only kinetic but also chemical equilibration starting from an off-equilibrium situation, using 
the transition-probability matrix elements of the underlying quantum field theory. 


2. Linear quark-meson a model 

To investigate the feasibility of a kinetic description of the off-equilibrium dynamics of the 
chiral phase transition the most simple two-flavor chiral model, based on the chiral group 
SU(2 )l x SU(2)r is considered, using a flavor doublet of Dirac fields ip, describing u and d 
quarks and antiquarks and a four-dimensional real-valued set of scalar fields (cr, n) transforming 
under the SO(4) representation of the chiral group 


15 . The Lagrangian reads 


C = ip[\$ - g{cr + iys p- r)]ip + ^(d^ad^cr + ■ d^n) - U(a,n) 


( 1 ) 


with the meson potential 


U(a,Tf) = ^-(cr 2 + vf 2 - u 2 ) 2 - Umlcr - U 0 , (2) 

where g 6 [3.3, 5.5] denotes the Yukawa coupling between quarks and mesons, A 2 = 20 the meson 
coupling constant (corresponding to a o mass of m„ ~ 600 MeV), f n the pion-decay constant, 
and v 2 = f 2 — m 2 /X 2 . The potential ([ 2 ]) contains the explicit breaking of the chiral symmetry 
due to the finite current quark masses resulting in a non-zero pion mass, m n ~ 138 MeV, of the 
pseudo-Goldstone modes tt. The constituent quark masses are given by m 2 = g 2 a q. 

The grand-canonical potential in mean-field approximation reads 

H(T, g) = U(a,n) + n^ (3) 


= ~ dn J (^3 \- E + Tln ( 1 + exp(-/3(£ - g))) + Tln(l + exp(-/3 (E + g )))], (4) 


with 








Figure 1. The phase diagram for the linear a model in mean-field approximation at /tb = 0, 
( 7 r) = 0. Left: the order parameter (a) and the effective cr mass m a = <9 2 f2/<9cr 2 . 


where E = p 2 + g 2 (a 2 + 7? 2 ) and the quark-degeneracy factor d n = 2NfN c = 12. The mean 
fields have to be evaluated self-consistently from the equilibrium condition 


<912 <912 
da <92? 


( 5 ) 


In the following we restrict ourselves to vanishing pion mean fields. A nice feature of this model is 
that by varying the Yukawa coupling g, one finds different kinds of phase transition as illustrated 
in Fig. |TJ 


3. Semiclassical particle-field dynamics 

The challenge in applying the above model to an off-equilibrium dynamical simulation of a 
system of particles (here quarks and antiquarks) and mean fields (representing the mesons) is 
that in order to reproduce the equilibrium-phase structure as depicted in [T] as the stationary 
limit, one has to ensure that both kinetic and chemical equilibration is possible through the 
introduction of the appropriate elastic collision terms for qq and qq scattering as well as quark- 
number changing processes such as qq -H- cr. In a full kinetic approach this is achieved by a set 
of coupled Boltzmann-Vlasov equations, which read in our case schematically (again restricting 
ourselves to the case of vanishing pion-mean fields) 

□cr + A(<7 2 - u 2 )a - f-n-m 2 + g (ipip) 

d* + -Jr ■ Vf- VgE^(t,x,p) ■ Vp-) 

Here, I and C denote collision integrals contributing to the meson-mean-field and quark-phase- 
space distribution functions respectively. 

In the following a novel scheme to Monte-Carlo simulate such a system of kinetic equations 
is defined, where one describes the mesons solely with a mean field and the quarks and 
antiquarks in terms of test particles. While the elastic-collision term is realized in a straight¬ 
forward way using the corresponding cross section from the underlying linear a model, 
one has to find a way to realize the interactions qq -H- a in such a scheme, while still 
fulfilling energy-momentum conservation and the principle of detailed balance, which are the 
fundamental principles constraining the off-equilibrium dynamics and ensuring the proper 
(Maxwell-Boltzmann) equilibrium limit. 

In our recently developed model (Dynamical Simulation of a Linear Sigma Model, DSLAM) 
this challenge is solved as follows: In order to properly simulate the collision terms on the right- 
hand sides of Eqs. 0 and 0 we define a space-time grid. In each time step at each spatial 


= I(a<^qq), (6) 

f q (t,x,p) = ->• <r>qq). (7) 
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Figure 2. While the energies of the 
quarks and the meson field in the thermal- 
box simulation show anti-correlated thermal 
fluctuations of the order A E/E 10 3 for 
the quarks and A E/E 10 2 for the field 
the numerical fluctuations of the total energy 
amount to only A E/e < 5 • ICG 5 . 


Figure 3. The total quark number in 
the thermal-box simulation shows fluctuations 
due to the dynamical pair-creation and 
annihilation processes, leading to an exchange 
of energy between the quarks and the mean 
field as shown in the left Fig. [2j 


cell the cross section for the annihilation process qq —>• a is used to stochastically determine an 
energy-momentum transfer from the initial qq pair, located in the cell. This energy-momentum 
change is transferred to the a mean field in terms of an appropriate relativistic Gaussian wave 
packet 5cr(t,x) (which simulates the gain term in /), and the qq pair is taken out of the test- 
particle ensemble (which simulates the corresponding loss term in C ). To simulate also the 
appropriate decay process, a —»• qq, we have to “particlize” the mean field locally in each spatial 
cell. This is done in the spirit of a course-graining procedure: First, the total energy-momentum 
content of the a field within the cell in terms of the corresponding cr-field energy-momentum 
tensor is determined. Then one assumes a local thermal equilibrium phase-space distribution, 
equivalent to this energy-momentum tensor. In order to fulfill detailed balance, the temperature 
has to be the same as that for the corresponding procedure for the quarks and antiquarks. The 
temperature is related to the mean-field value which depends on the scalar quark-antiquark 
density. In this way a temperature can be determined. It is important to note that it is defined 
in the local rest frame of the heat bath and thus the cr-phase-space distribution is given by a 
Maxwell-Jiittner distribution f a oc exp[— p^u^{t, x)/T], where u M is the four-velocity, given by 
the total field-four-momentum in the spacial cell under consideration, u^ = p^/E. Now in each 
time-step within each spatial cell one can choose an ensemble of a particles according to this 
local Maxwell-Jiittner distribution and using the corresponding qq —> a decay rate to determine 
the gain term to C. The loss term for the mean-field equation in the collision term / is again 
achieved by taking the appropriate amount of energy and momentum out of the mean field in 
terms of a Gaussian wave packet. 

In summary we have achieved a scheme which enables us to simulate the set of Boltzmann- 
Vlasov Eqs. © and 0 using test particles for the quarks and antiquarks and restricting the 
description of the mesons strictly to the mean-field level. The scheme by construction fulfills 
energy-momentum conservation through the Gaussian wave-packet description for the exchange 
of energy and momentum between the mean field and the test particles. At the same time 
also the principle of detailed balance is fulfilled, using the coarse-graining approach to locally 
map the field-energy-momentum distribution to a local-equilibrium Maxwell-Jiittner distribution 
to reinterpret the mean field as a meson phase-space distribution and using the leading-order 
transition rates (cross sections) of the underlying QFT linear a model fulfilling the detailed- 
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Figure 4. The distribution of the quarks’ 
kinetic energy in a thermal-box calculation 
show an excellent agreement with the expected 
Maxwell-Boltzmann distribution. 



Figure 5. The distribution of the kinetic field 
energy to its Fourier modes. The classical “UV 
catastrophe”, i.e., the equipartition theorem 
associating an average energy of k^T/2 for 
each Fourier mode, is avoided by the use 
of a finite finite width a e g of the Gaussian 
wave packets for the transfer of energy and 
momentum between particles and field. 


balance principle. 


4. Proof of principle: “Box” calculations 

As a plausibility check of the simulation method the stability of an equilibrium situation has 
been tested in a finite cubic box with periodic boundary conditions 14 . The stability of the 
energy conservation is demonstrated in Fig. [2j While the kinetic energy of the quarks and the 
energy of the mean field show anti-correlated thermal fluctuations the total energy stays stable 
within the numerical accuracy of the simulation. 

The annihilation and creation processes of quark-antiquark pairs lead to thermal fluctuations 
in the total particle number, as shown in Fig. [3j 

The energy distribution of the quarks shows the expected Maxwell-Boltzmann distribution 
at the expected temperature, as demonstrated in Fig. [4} Also the spectral analysis of the mean 
held shows that the thermal equipartition theorem for the kinetic held energy is fulfilled at 
small wave numbers (long wavelengths), i.e., each mode contains an average energy of k^T/2. 
On the other hand the “UV catastrophe” must be avoided due to the finite total-energy content 
within a finite box. Indeed, the energy-momentum transfer between particles and the held due 
to the pair-annihilation and -creation processes is not strictly “local” but occur within a hnite 
volume whose scale is hxed by the hnite width, a, of the Gaussian wave packets used as held 
increments to keep care of the correct energy-momentum transfer in each process. Thus the 
large wave numbers (short wave lengths) are effectively cut off at a scale fc cuto ff = V^/cr, as 
shown by plotting the corresponding Gaussian on top of the /c-distribution of the kinetic held 
energy (Fig. [5]). 

The dynamical behavior of the held in our thermal-box simulation is illustrated in Fig. |hj 
showing the held distribution within the xy plane of the simulation. Starting the simulation 
with a uniform mean-held equilibrium value, after a short time some local blob-like disturbances 
have developed (left panel) due to quark-pair-annihilation processes, leading to the propagation 
of Gaussian wave packets on top of the still quite uniform average held value, as implemented by 
our concept of describing the exchange between held and particles in these processes. The right 
panel shows the fully developed equilibrium-thermal-held fluctuations in the long-time limit of 

































Figure 6. Starting with a uniform 
mean field, a short time after the start 
of the simulation the held is disturbed 
by some Gaussian wave packets due to 
quark-pair annihilation processes (left 
panel). In the long-time limit the full 
equilibrium-thermal-held fluctuations 
have developed. 
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Figure 7. The held fluctuations are 

Gaussian distributed around an average value 
which can slowly drift with time due to the 
dynamical fluctuations of the held momentum. 


the simulation. In this limit the held values are Gaussian distributed around a mean value (cf. 
Fig-0, which is expected due to many random energy-momentum transfers between the held 
and particles. Thereby the average held value can slowly drift with time due to the thermal 
fluctuations of the held energy and momentum. 

5. Conclusions and outlook 

In this talk it was demonstrated that the dynamical description of the chiral phase transition 
in a simple quark-meson model is feasible with a novel Monte-Carlo-simulation technique for a 
corresponding coupled Boltzmann-transport equation for the meson mean-held and the quark 
and anti-quark phase-space distribution functions, implementing both elastic quark and anti¬ 
quark scattering as well as quark-antiquark-pair creation and annihilation processes, enabling 
both kinetic and chemical equilibration between particles and the mean held (mesons). 

The simulation is set up in a way that both energy-momentum conservation and the principle 
of detailed balance are precisely realized (within the limits of achievable numerical accuracy). 
While the elastic quark-scattering processes are simulated with a straight-forward test-particle 
realization, the particle-held kinetics has demanded the development of a novel scheme. 

The qq ^ cr annihilation process is evaluated by the Monte-Carlo sampling according to 
the corresponding transition matrix elements from the underlying quantum-held theoretical 
interpretation of the linear a model. The corresponding energy and momentum are precisely 
transferred to the a-mean held in terms of an appropriate disturbance in form of a relativistic 
Gaussian wave packet. To obey the principle of detailed balance, also the inverse a —> qq decay 
has to be simulated. To that purpose a coarse-graining procedure based on the local energy- 
momentum content of the held has been used for a “particlization” of the mean held in terms of a 
local Boltzmann-Jiittner equilibrium distribution, which in turn enables a Monte-Carlo sampling 
of the cr-meson decay according to the decay rate from the quantum-held theory picture of the 
model, which automatically implements detailed balance. 

This scheme has now been applied to off-equilibrium situations as in a “thermal quench” 
in a box where the particles and helds are initialized with different temperatures and it was 
demonstrated that the distribution can describe a possible phase transition from the initial 
state to a late-time (equilibrium) state. 

Last but not least also the case of “expanding fireballs”, mimicking the situation of the 
medium created in heavy-ion collisions, is under study. 
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